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SUMMARY 


The general equations of motion for an airplane with a number of 
spherical fuel tanks are derived. The motion of the fuel is approxi- 
mated by the motion of solid penduliams. The same type of derivation 
and equations are shown to apply to any type of fuel tank where the 
motion of the fuel may be represented in terms of undamped harmonic 
oscillators. 

Motions are calculated for a present-day high-speed airplane and a 
free-flying airplane model with two spherical tanks in the symmetry plane. 
These calculations show that the normal airplane motion may be consider- 
ably modified and that residual oscillations may result. The ratio of 
the natural fuel frequency to the natural airplane frequency is shown to 
be the most important parameter for determining the effect of the fuel 
motion on the airplane motion. The stabilizing effect of turbulence in 
the fuel is discussed, and it is suggested that the stabilizing effect 
of artificially induced turbulence be investigated experimentally. 


INTRODUCTION 


Small-^plitude lightly damped lateral oscillations are a trouble- 
some characteristic of certain high-speed airplanes. Several possible 
explanations for these oscillations, which are adequate in specific 
cases, have been offered. For example, reference 1 shows that nonlinear 
aerodynamic derivatives could cause such oscillations, and it has been 
shown that atmospheric turbulence is another possible cause. It has also 
been suggested that a possible cause of such oscillations is the motion 
of fuel in the tanks. In some recently designed airplanes the mass of 
the fuel relative to the airplane mass is much larger than has been 
common in the past; therefore, the effects of fuel motion can be expected 
to be relatively more important. In fact, in several cases baffling the 
fuel tanks was found to have considerable effect on the general handling 
qioalities of the airplane and sometimes actually eliminated the trouble- 
some lightly damped lateral oscillations which had been present. 
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An experimental investigation of the effects of fuel motion on the 
lateral motion of a free-flying airplane model is described in refer-;' 
ence 2 . The results indicated that the effects of fuel motion were 
noticea'ble and caused the lateral motion of the model to he very erratic. 

The present analysis treats each fuel tank as a pendulum oscil- ' 
lating in two degrees of freedom and applies Lagrange's equations of 
motion to obtain the interaction between these pendulums and the air- 
plane. Thus, for small motions the fuels are treated as simple harmonic 
oscillators. The results are applied to obtain the general equations of 
motion of this system and, in particular, the lateral motion of an air- 
plane with internal fuel tanks in the plane of symmetry of the airplane. 
Since the general solution of the equations is extremely complicated, 
an attempt is made to evaluate the results by carrying out numerical 
calculations for specific cases. This approach is shown to be adequate 
in yielding the most general effects of fuel motion. 

The discussion of the numerical application of the equations of 
motion to specific cases is given in detail following the derivation of 
the equations of motion. This discussion of results is understandable 
quite independently of the derivation of the equations of motion. 


SYMBOLS 


X, y, Z airplane stability axes with origin determined by 

equations (l 3 )j also components of applied forces 
along these axes 

L, M, N components of applied moments about X-, Y-, and 

Z-axes, respectively 

unit vectors along X-, Y-, and Z-axes, respectively 

components of translational displacement of airplane 

vector translational velocity of airplane 
( ix + jy + kz = + v) 

Uo magnitude of steady-state velocity 

w components of disturbance translational velocity of 

airplane 

vector disturbance velocity of airplane 
(_iu + j^v + kw) 


h A> ^ 

X, y, z 

* 

r 


V 
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R 


V 

3 

0, e/i 


cu 


T-, 


I 

h 

k 


, vector position of a point in airplane 
(l?x + + iE^z) 

vector position of 'center of gravity of fuel in 
a particular fuel tank 

. total vector velocity of a point in airplane 
sideslip angle 


(tan~^ — ^ 

V ^oj 


inf initesmal rotations of airplane about X-, Y-, 
and Z-axes, respectively 

vector rotational velocity of airplane 

( i0 + 10 + 


components of angular displacement from vertical of 
line joining fuel center of gravity -to tank center, 
taken on mutually perpendicular planes; ^ positive 
in direction of positive roll and t) positive in 
direction of positive pitch 


distance from tank center to fuel center of gravity 

vertical displacement of fuel center of gravity from 
equilibri\jm position 

number of fuel tanks 


m 


mass 


®t 


total mass of airplane and fuel 




} 


total moments and product of inertia of airplane 
about X-, Y-, and Z-axes 



rigid-body moments and products of inertia about axes 
through center of gravity 


fuel moments of inertia about and T)-axe s' through 
tank center 





k 
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If = 


% 


Kz 


%Z 


Kf = 


I Z 


Ek 

®P 

P 

Il/2^ I'2 


t ' 
s 


g 
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P 

S 


fuel moments of inertia atout and T]-axes through 
fuel center of gravity 


nondimens ional radius of gyration in roll 



nondimens ional radius of gyration in yaw 



nondimens ional product-of-inertia parameter 



kinetic energy 
potential energy 
period 

time for exponentially damped or increasing oscilla- 
tion to. halve or double amplitude, respectively 

time 

nondimens ional time parameter 
acceleration due to gravity 

flight-path angle with respect to horizontal 
air density 
wing area 


¥ 
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b 

*^bf 


Cl 


wing span 

lateral nondimensional mass coefficient 

ull \ 

differentiation operator 

cos 7 o> 


trim lift coefficient 


I pUo^S 


rolling-moment coefficient /Rolling moment^ 


I pUo^Sb 


"n 


yawing -monient coefficient 


^Yawing moment^ 

I pUo^Sb 


Cy 



lateral-force coefficient 


/Lateral force^ 


I 


I PUq^S 
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particular fuel tank, or summation index over fuel 
tanks (f = 1, 2, . . . k) 

airplane without fuel 

initial conditions at t = 0 


DERIVATION OF EQUATIONS OF MOTION 
Assumptions for Derivation of Gteneral Equations of Motion 


As a first approximation, only the effect of the motion of the fuel 
as a whole is considered; that is, only the fundamental mode of the wave 
motion is considered, and this mode is approximated by rigid-body 
motion. The main effect of the internal wave motion is to introduce 
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damping inio the fuel oscillation. This damping is caused by the con- 
version of kinetic energy into heat through the txirbulence caused by 
the splashing of the fuel. A strictly analytic consideration of such 
damping effects is extremely difficult; on the other hand, the damping 
caused by the viscous tangential forces between the fuel and the tank 
is completely negligible (see reference 3 ). The analysis of the problem 
is therefore confined to the motion with no fuel damping and the effect 
of the damping is considered in the discussion of the results. 

In a spherical tank the fuel can oscillate approximately as a rigid 
body if no splashing is assumed for- small oscillations. The motion may 
be pictured as the "rocking" of a spherical segment of constant shape. 
The restraining force of the tank, which always acts in a direction 
normal to the motion, is exactly analogous to the tension in a pen- 
dul\mi. Thus, the small motions of the fuel in a spherical tank may be 
represented quite well by the well-known. simple properties of small 
pendulum motions. This approach is used in the mathematical analysis 
of the problem. 

The effect of aspherical tank shape can be approximated by 
replacing the tank by an equivalent harmonic oscillator with an arbi- 
trary amount of turbulence damping added even for small motions. For 
example, the representation of rectangular tanks- as harmonic oscillators 
is discussed in reference 4 . Thus in this case also the most general 
effects of the fuel motion on the airplane motion should be qualita- 
tively obtainable by this type of analysis. 

The effects of large -amplitude fuel motions will be discussed 
qualitatively after the discussion of the results of the mathematical 
analysis. As usual in stability analysis all motions are assumed small 
and second-order terms are ignored. 


Derivation of General Eqioations of Motion 

With the preceding assumptions the physical problem can be con- 
sidered as the interaction between two or more rigid bodies, namely the 
airplane and the several fuel pendulums, with each fuel pendulum con- 
sidered as suspended from the tank center. The only potential energy 
considered in the system is that of the pendulums. ^ if the inertial 
characteristics of the airplane and the fuel are known, the kinetic 
energy of the system can be obtained from the translational and rota- 
tional velocities of the airplane and the fuels. With this information 
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the interactions in the system can be obtained by using Lagrange’s 
equations of motion in the form (see reference 5) 







8qi 8q^ 


Qi 


(i = 1, 2, 


( 1 ) 


n) 


where qj^ is one of the n generalized coordinates of the system, cor- 
responding to the n degrees of freedom, qj^ is the corresponding 
velocity, and is the corresponding generalized force. The qj^ 

will be lengths and angles and the corresponding will be forces 

and moments, respectively. 


The airplane Itself introduces the customary six degrees of free- 
dom, which are the three displacements of the airplane system along 
axes fixed in the airplane (x, y, z) and the corresponding angles of 
rotation of the airplane about these axes (0, 0, ilr). For small dis- 
placements the pendulim motion can be described by two angles, ^ 
and T|, since the vertical motion can be neglected (see fig. l). The 
angle ^ ' is measured from a vertical line through the tank center to 
the projection of the line joining the tank center to the fuel center 
of gravity on the vertical plane parallel to the Y-axis and q is the 
corresponding angle in the vertical plane parallel to the X-axis. For 
small angles, ^ and q may be represented as in figure 1. In effect 
this figure makes use of the fact that small angles may be added vec- 
torially. When the two additional coordinates ^ and T| are used to 
describe the pendulum motion, the whole system has two additional 
degrees of freedom for each fuel tank. 


Expressions must be obtained for E^ and Ep in terms of the 
coordinates of the system and their time derivatives in order to use 
equation (l). The only potential energy is that of the fuel pendulums', 
which can be written as follows for k fuel tanks 

k 

®P = IZ “f^S (2) 

f=l 


For the height of the center of gravity, in each fuel tank, as can be 
seen from figure- 1, 


h = 


+ h^ Z - 2 cos ^ cos q 
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or . 

h«.|z(^2 + ^2) (3) 

Note that the vertical displacement h is of second order in the small 
quantities T| and This fact* justifies the previous statement that 

the vertical displacement cotild be neglected in describing the pendulum 
motion only by the two coordinates tj and As might be expected, 

equations (2) and (3) indicate that each fuel pendulum is being con- 
sidered as an undamped oscillator with two degrees of freedom in a 
horizontal plane. 

The kinetic energy of the total system can be. written as the sum 
of the kinetic energy of the airplane and the kinetic energies of the 
fuels. Also the kinetic energy of each rigid body can be expressed as 
the sum of the translational energy of the mass moving with the velocity 
of its center of gravity and the rotational kinetic energy of the mass 
about its center of gravity. Thus, when the inertial characteristics 
of the airplane and the fuels are known, the kinetic energy can be 
obtained as a function of' the generalized coordinates and velocities if 
the translational velocity of each center of gravity and the angular . 
velocities of the airplane and fuels about their respective centers of 
gravity can be expressed in terms of these generalized coordinates and 
velocities. 

In order to obtain the required expressions for these velocities 
a system of axes fixed in the airplane with the X-axis along the steady- 
state velocity at t = 0 is used, as is customary in stability analysis. 
For the present the origin of the coordinates will not be specified. 
However, these stability axes are not inertial axes and Newton’s second 
law applies only in an inert ia .1 system of axes. The inertial axes may 
be taken as axes fixed in the earth. Then in the equations of motion 
the velocities and accelerations must be measured with respect to the 
earth, and their expressions in terms of components in the moving air- 
plane axes may be obtained as shown in reference 6. These expressions 
will -give the kinetic -reaction forces, which for the case of a rotating 
system are often referred to as "gyroscopic" forces. For the velocity, 
referred to the inertial system, of any point . defined by the vector R 
in the airplane axes ( in particular, for the centers of gravity pre- ~ 
viously discussed) 


y = r + a)XR + R= ^Uq + v-f^xR + R 


( 4 ) 
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where all vectors are given in terms of the airplane axes and 

r = + kz = ^Uq + v 


is the velocity of the origin of the airplane system with respect to 
the earth, while v and o) are the translational and rotational dis- 
turbance velocities of the airplane axes. 

Equation (4) may now be used to express the inertial velocities of 
the airplane center of gravity and the fuel centers of gravity in terms 
of the generalized coordinates by inserting for R the values Rq 
and R„, where Rg^ is the vector position of the airplane center of 

- I 

gravity and R^ indicates the vector position of any particular fuel 
center of gravity. The vector ^ is constant; therefore, ^ = 0. 

To obtain R^, note that to first order 


5f «« Sf ■ ® " ’Ifl ■ J.^f(^f ■ ^ cos,7o “ ^ sin 7 q) + 

klf cos^7o + 0 - 

5 f ss - iif[sin 7o + (e - 'nf)cos 7p] - cos 7 q - t sin 7 q) + 

kZfjcos 7 q - - qjjsin 7 qJ 

where R^ is the fixed position of the tank center. Since 7 q is 
constant, to first order 

” -i^f (^f - ^ cos 7 q - ^ sin 7^) + 

kZf - e)sin 7o 
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Again keeping only first-order terms leads to the following eqiiation: 

^XRfRja)X]^ + cos 7o - j.^f cos 7o + t sin 7o) + sin Jq 

Now combining the last two equations gives 

^ X ^ X + iq^Z^ cos 7o ~ + kZ^q^ sin Jq (5) 


This equation shows, as could be predicted physically since no viscous 
force is assumed between the tank arid the fuel, that the airplane rota- 
tion affects only the motion of the tank center. 

From equations (4) and (5) the necessary translational velocities 
can be obtained for the translational kinetic energies. The rotational 
velocity of the airplane is simply oj. The spinning motion of the fuel 
about the vertical axis is ignored;, then, the rotational velocity of 
the fuel may be given by the componeiits ^ and q. The two corre- 
sponding horizontal axes ‘of rotation through the fuel center of gravity 
are principal axes of the spherical segment of fuel; consequently, no 
product -of- Inertia terms occur in the fuel rotational energy. Also, 
since the airplane center of gravity is in the airplane symmetry plane, 
lyz " ^XY ” ® only the product of inertia will appear in 

the airplane rotational energy. 

By use of equation (4), the airplane velocity can be shown to be 
Ya = i(Uo + u + 0Rz^ - tRya) + j(v + tRxa -,^Rza) + 


k(w + 0Ry^ - 6RxJ (6a) 

When eqiiation (5) is substituted in equation (4), the velocity of any 
particular fuel center of gravity is 

Vf = i(Uo + u + RR^^ - +■ qf Zf cos 7o) + + tRx^ - Zf] + 


k(w + 0Ry^ - RRx^ + flf-Zf sin 7o) (6b) 
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If V is the magnitude of the translational velocity of the center 
of gravity of a rigid body, l' its moment -of -inertia tensor about the 
axes through its center of gravity, and w the rotational velocity of 
the rigid body, the kinetic energy is 


= i mV^ + 

K 2 




X* 2 x> 2 


2I^Z“y“z - 2Ixz“x‘^z " 


( 7 ) 

Thus, for the kinetic energy of the airplane, substitution of equa- 
tion (6a) in equation (7) gives 





+ u + 6Ry 






ya 


- 6R 


'^f~\ 






-2 

2 ^ 




and, for the kinetic energy of each fuel, substitution of equation (6b) 
in equation (7) gives 




+ u + 0R„ 

Zf 


+ T|fZf cos 7o)^ + 


^v + ij'Rx^ - ^Rzf - ^f ^fj^ + + ^Ry^ - SRxf + hf Zf sin 70 






(8b) 


k 

For the total kinetic energy, E]^ = Ej^^ + y therefore, 

f=l 

equations (2), (3)^ and (8) may be used in equation (l) to obtain the 
(2k + 6) equations of motion. However, it must again be recalled that 
the coordinate system is rotating. The whole system is therefore sub- 
ject to an additional gyroscopic acceleration since the time-derivative 
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operator contains an additional gyroscopic term (see reference 6) when 
the components of the velocity (or any vector) are taken in the rotating 
system; 



r + CD X r 


Thus the gyroscopic 
CD X r = 


acceleration acting on the whole rotating system is 
03 X |iUo + v) (^ X i)Uo'= JijfUo - k6Uo 


The effect of this acceleration can he brought into the equations of 
motion by considering the inertial reaction of the total mass to this 
acceleration as an additional applied force. If the total mass is 
k 

iiLj. = ma + y mp, this reaction has the following components: 
f=l 


Y- = - mpUot 

Z* = mpU^e 



(9) 


In addition there is the inertial reaction torque ^ on the fuel; this 

torque acts about the tank center and is caused by the acceleration of 
the tank center. For each tank, the vector reaction torque is 


Mf ~ (Sf “ %) ^ ™fUo(“i''^ 


= nipUQ y^-iZp sin 7 q + k^f cos 7^) x + kfi ) 

sa nifZfUo^iiir cos 7 q + ^6 sin 7 q + kij^ sin 
Since m’ = cos 7q + sin 7q and = 


+ Second-order terms 


n 


j 
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li<- 


M 


^ = mfZfUo\if(cos^7Q + sin^ro) = mfZfUoiif 


(10a) 


= nifZfUoe sin 7 q 


(10b) 


The forces and moments in equations (9) and (lO) must be added to 
the weight and aerodynamic forces to obtain the Qj^ in equation (l). 

For convenience, the results of equations (2), (3)^ and (8) are 
as follows: 


Ep - i g y mf Z^. 


(11a) 


f=l 


Ek = ^[(Uq + - iRy^)2 + (v + tRx^ - ^RzJ2 + 


k 

(w + ^Ry^ - ORxa)^ + XI 
' J f=l 


nij^ 

~2 


[(Uo + u + 6Rzf - i^fRyf + Tlf Z^ cos 7 q)^ 


(» + - ((E^^ - ifhY * (“ + + EfJf sin r„y 




Eq\iations (9), (lO), and (ll) may now be used in equation (l) to 
obtain the equations of motion. It should be noted that ' 

— - — 7 — ^ ^ = — , — = — , and — = — . For example, to obtain 

8x 8(Uo + u) 8u Sy Sv 8z 8w 

the equation of motion in the x-direction^ note that 

^Ek _ _ 0 

Sx. dx 
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Then, equation (l) may be written as follows: 


dt \^x / dt / 


= ma(u + 6Rz^ - ^ (li + ^Rzf " iRyf + V,^ 2^ cos 7o) 

= m^U ' + e^m^Rz^ + ^ m^Ry^^ + 

k 

.cos 7o ' ( 12 ) 

The position of the origin of airplane coordinates has not yet been 
specified. Equation (12) and the similar eqioations obtained in the 
other degrees of freedom suggest that the position of the origin be 
determined by the following three conditions: 


“a«x_ + Y1 Rxf = 0 

f=l ^ 

“aSya ^ ^ ' ° 

“a"z, + ^ mj.E. . 0 
a f=l ^ 


( 13 a) 


(13b)- 


(13c) 


Equations (13) imply that the origin is at the position of the total 
center of gravity when the fuel mass is treated as being concentrated 
at the tank center.. This choice of the origin greatly simplifies 
expressions such as equation (12). The physical reason for this choice 
is again that the fuel does not rotate with the airplane; thus, a force 
acting on a line through this point, the center of gravity where the fuel 
reaction is assumed concentrated at the tank center, will produce no 
rotation of the airplane. 
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The following substitutions will also greatly simplify the writing 
of the final equations of motion: 


% » 4 * ”a[(%a)" " ^ [(«7f )" * 

Iz = Iz * “a[(Bxa)^ * (“ya)^J * ^ 

X 

^XZ =. ^XZ ^ 


> (l^a) 


Note that the quantities defined by equations (l^) are the total moments 
and product of inertia about the origin of the airplane coordinates 
when the fuel mass is assumed to be concentrated at the tank center. 
Finally, the necessary moments of inertia of each fuel pendulum about 
the tank center are 




Without loss of generality then, equations (13) and (l 4 ) are used in 
the equations of motion obtained by substituting equations (9), (lO), 
and (11) in equation (l). The general equations of motion can now be 
given as follows: 




+ mfZf2 


(l 4 b) 
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m+ii + cos 7 ^ = X 

f=l 

jj^ 

®t(* " ^ lOfljTif Sin 7o = Z 


Iy0 - - ^%Rz^Ry^ + nifR^^y^il 

^ (Rzf cos 7o - Rxf sin 7o)“>f “ 


+ u cos 7o + (w - Uoe)sin 7 q + sin 7©. + 

S^Rzj^ cos 7 q - Rxj^ sin 7©) - Ry^'*' ° 
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In equations (l5a) and (l5t>) only the fuel equations for the first fuel 
tank have been written. In each set there are k similar fuel equa- 
tions. As has been previously stated, the forces on the right-hand sides 
of these equations are the applied forces and the weight and aerodynamic 
forces. 


Simplifying Assumptions 

The equations of motion have been separated into what would gen- 
erally be considered the longitudinal motions, equations ( 15a), and 
the lateral motions, equations (l5b). In the ordinary six-degree-of- 
freedom case, as can be shown frOm considerations of symmetry, no 
cross-coupling terms exist between these motions in the aerodynamic 
forces (see reference 7). Although such terms are known to exist in 
practice, they are small and generally neglected. However, many cross- 
coupling terms occur between equations (l5a) and equations (l5b) because 
of the fuel motions, even when the aerodynamic coupling forces are 
ignored. The magnitudes of these fuel forces can be seen to depend oh 
the masses of the fuels, the vector positions of the fuel tanks, the 
"pendulm length" (i.e., the radius of the tank and the height of the 
fuel), and the accelerations involved in the fuel motion. Increasing 
the magnitude of any of these parameters will increase the effect of 
fuel motion. For most airplanes these fuel forces will be relatively, 
small, but the present investigation is primarily concerned with all 
first-order fuel effects. In any particular case, the actual magnitudes 
of these forces may of course be obtained by inserting the values of 
the previously mentioned parameters. 

Symmetrical fuel distribution .- Some of the terms in equations (l5a), 

for example, the 0 and ijr terms in the 6 equation, are essentially 
product-of-inertia terms arising from unsymmetrical distribution of the 
fuel about the y = 0 plane (i.e.-, the plane of symmetry). That this 
is so can be seen if the fuel tanks are assvuned to be distributed 
symmetrically with respect to the symmetry plane; that is 


k 

r “fHxjRyj 


k 

= r 

f=i 




k 

n “"f%f = 0 

f=i ^ 


(I6a) 


and from equations (I 3 ) 



(16b) 
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therefore, these terms vanish for symmetrical fuel distributions. In 
most cases the fuel will be symmetrically distributed, and substitution 
of equations (l6) in equations (15) yields the equations of motion for 
symmetrical fuel distribution: 


m^u + ^ cos "/q = X 

f=l 


m 


- Uo^) + Y 1 7o = z 


ly® zZ. (^zf M 


I 

+ “l^l[s\ + u cos 7o + (w - Uo0)sin + Ry^‘^ sin y^ 
cos 7o - Rxi sin 7o) - Ryj^'t cos 7o] = 0 


> (iTa) 


/• • \ 

m^fv + UQ^t' ) - L = Y 

' . f=l 

Iz^ - lx# - ^ ”^^f (^^flf + RyfRf cos. 7o) = N 


K 

“ Ixz^ ^ |Z “f^f(^Zf^f + ^^iif sin 7q) = 1 


> (17b) 
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In eqxaations (17) even though the terms arising from unsymmetrical fuel 
distribution have vanished^ some cross -coupling terms still remain 
between equations (l 7 a) and (l 7 b). These terms occur in the t] equa- 
tions of set (l 7 a) and in the 0 and \|/ equations of set (l 7 b). The 
significance of these terms is evident since each contains a factor Ry^. 

Thus, these terms arise when the airplane has fuel tanks with centers not 
in the plane of symmetry, even though they are symmetrically distributed 
with respect to this plane. For example, they would arise for wing-tip 
tanks. Physically, these terms clearly give the interaction between' 
the longitudinal fuel motion T|f in the wing-tip tanks and the airplane 

rotation about the vertical axis, which consists of the lateral 
motions 0 and t* For example, assume for simplicity that 7 q = 0 ; 
then, a yawing acceleration of the airplane will cause a longitudinal 
fuel acceleration *1]^ in the wing-tip tanks, and vice versa. 

From this discussion the motion appears to couple the lateral 
and longitudinal airplane motions even for the perfectly symmetrical 
fuel distribution described by equation (l6a). However, the fact that 
this coupling does not occur can be seen by considering any pair of 
symmetrically placed and loaded tanks. Designate the. q motion in 
this pair of tanks by and q^. Then, the q eqviations in equa- 

tions (l7a) show that a longitudinal horizontal acceleration of the 
tank gives rise as expected to q accelerations. Since the system 
is linear this portion of the q motion may be considered independently, 

and because of the symmetry of the two tanks it' is seen that Ti-, = 'ti 

*1 *2 

for the portion of the q motion arising from the longitudinal motion. 
Therefore, ■ in the \tr and 0 equations of set ( 17 b) the effects of 
this q will vanish since Ry^ = .-Ryg. In a .similar manner the laterally 

caused q motion can be shown to have no effect on the longitudinal 
motion. Essentially the argument is that the q motion for each pair 
of tanks can be split up for perfectly symmetrical fuel distribution 
into symmetrical and antisymmetrical motions. The symmetrical portion 
of the q motion for each pair of tanks couples only with the longi- 
tudinal motion; the antisymmetrical q motion couples only with the 
lateral motion. Thus, for perfectly symmetrical fuel distributions 
the q equations of set (l7a) could be combined with, set (l7b), 
only lateral degrees of freedom in the q equations being used ( since 
the symmetrical portion of the q motion .is of no interest); or 
the q equations may be used as shown with equations (l7a), the lateral 
degrees of freedom in the q equations being ignored. 

Fuel tanks centered in symmetry plane . - Many airplanes have' large 
internal fuel tanks which are centered 'in the airplane symmetry plane. 

For such airplanes the equations of motion may always be separated into 

independent lateral and longitudinal modes since R^ = 0. , Using this 

^f 
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value in equations (17) gives the equations of motion for tanks- centered 
in the symmetry plane: 


mtii + ^ COS 

f=l 




“tl* ■ ^O®) + II sin 7o = Z 

f=l 

lY^- . t K COS 7q - Rx^ sin = M 

+ “lh[sni + u cos 7o + (w - Uo6)sin 7 q + 
cos 7o - sin = 0 

# 

mt(v + Uoir) - ^ = Y 

f=l 

.. k 

^ 

f=l ^ 

k 

- hz^ + C %^fRz/f = L 

- V - Uq 4' - R^^’ii; - = 0 


> (18a) 


> (l8h) 
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Each set of equations, (l8a) and (l8b), contains (k + 3 ) variables; 
the two sets can be. seen to be independent of each other since the 
i motion of each tank couples with only the lateral motion and the'^ 

■q motion couples with only the longitudinal motion. 

For the case of a single fuel tank at the airplane center of .gravity 
the modification of equations (18) is obvious. Then " 5 f “ 

all coupling between the rotational motion and the fuel motion vanishes; 
that is, all the fuel terms in the rotational equations vanish and all 
rotational terms in the fuel equations vanish. For an aspherical tank 
the rotational coupling in this case will be small. 


Limitations Inherent in the Approximations 

Before proceeding to the application of equations (l8b) it is 
appropriate to consider somewhat more explicitly the assximptions involved 
in the indiscriminate dropping of all second-order terms which appeared 
during the derivation of the equations of motion. In this connection 
the correction, arising from. the airplane accelerations, to the constant 
acceleration fie'ld g Involved in the pendul'jm potential energy should 
be considered. The assumption which is implied in neglecting these 
accelerations is that the accelerations of the tank centers are small 
with respect to g. , 

Considering this and previous approximations, it can be seen that 
three essential assximptions were made in dropping second-order terms: ' 

.(1) The fuel and airplane angular displacement variables are small 
enough so the angle approximates its sine. However, this approximation 
sometimes took the form that the angle was much less than 1 radian. 

(2) The disturbance velocities are much less than Uq, and products 
of the linear or angular velocities can be ignored. 

(3) The accelerations of the tank centers must be small compared ' 
to g. 

Strictly speaking then,' the statement that the equations of 
motion (15), and also the simplified equations, are accurate eqiiations 
of motion to first order is to be taken to mean that the motions to 
which these equations apply are restricted by the preceding three con- 
ditions. Thus, the equations would appear to remain accurate at least 
at the beginning of a disturbance. Moreover, when the motion becomes 
large enough that these assumptions break down, the fundamental physical 
assumption that the fuel may be considered to move as a rigid body also 
breaks down; therefore, nothing can essentially be gained by keeping 
higher-order terms in the mathematical expressions. 
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Since the penduliom motion is little changed even up to angles 
of 30° to it could be expected that aside from splashing effects 

these equations should remain a good to fair approximation even at such 
angles. On the other hand, even the splashing effects, although they 
would introduce some damping and change the inertial characteristics 
of the pendulum somewhat, could certainly not be expected to cause the 
general assumptions to break down completely for fuel motions up to 
angles of 30° to 1*0°. Therefore, the equations of motion derived are 
assumed to present a fair picture of the disturbance motion even up to 
fuel displacements of this magnitude. 


APPLICATION TO SEVERAL CASES 
Nondimens ional Eqviations 'for Tanks in Symmetry Plane 


The equations of motion (l8b) have been applied to the lateral 
motion in several cases with two fuel tanks in the plane of symmetry. 

In these cases the lateral motion can be considered independently. The 
applied forces are the weight, the- usual aerodynamic forces linear in 
the disturbance velocities, and any distxirbing forces that may be 
present. In order to put the equations in nondimens ional form, the 
nondimens ional lateral airplane equations are used as obtained in 'refer- 
ence 8. The fuel equations are made nondimens ional by making the 
standard transformation to nondimens ional time derivative, as in the 

airplane equations, and then dividing through by — — . The resulting 

nondimens ional expressions in the following equations are defined in 
equations (20): 


- i Cyp^p + i Cl tan - 


^^b2 
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■ (^ D^ + d)^ + ^ D^0 + = 0 (I9d) 


-D3 - 1)2 + D^ ^If + 1)20 + + G^ ^ = 


0 ■ (I9e) 


where 


g.£% %=- 4 t 

Uo^ nif.Zpb 




( 20 ) 


The derivatives Cy and Cy were assumed to be zero. 

-^p -^r 


Methods of Solution 

In the present case the two fuel degrees of freedom introduce two 
additional oscillatory modes into the characteristic solution^ in 
addition to modifying the original airplane mode. The motion will 
therefore be a combination of three oscillations (aside from the leqs 
important exponential modes) ^ but just knowing the three oscillatory 
roots is insufficient to indicate the type of motion since the relative 
magnitudes of the oscillatory modes must also be^known. For this 
reason motions must be calculated in order to see the actual effects 
of the fuel motion. However^ in several cases the characteristic roots 
were also found in order to facilitate the interpretation of the motions 
These cases will be discussed subsequently. 
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The most convenient method for the analytical solution of a set • 
of linear ordinary differential equations such as equations ( 19 ) is 
probably the Laplace transform method (see reference 8 ). However this 
method is extremely cmnbersome and difficult to check since it involves 
the expansion of fifth-order determinants in which the elements are 
often quadratic functions of the characteristic root. Therefore, it 
seemed preferable to use some step-by-step method which would be more 
amenable to machine computation. 

Reference 9 gives a matrix method for getting the step-by-step 
solution of a set of linear ordinary differential equations. When 
applied in the present case to equations ( 19 ) this method results in 
a simultaneous solution for the motion in each of the five degrees of 
freedom and also for the motion in D0, Di|f, and The cal- 

culations were carried out on the Bell Telephone Laboratories X-66'jkk 
relay computer in use at the Langley Laboratory. The essential details 
of the method are given in appendix A. 


Solutions for Several Cases 

The two basic cases for which motions were calculated were case A, 
a present-day high-speed airplane with two fuel tanks satisfying the 
conditions .for equations (19)^ and case B, which corresponds essentially 
to case B of the model used in reference 2. The essential parameters 
for these two cases are given in tables I and II. Table II gives the 
conditions for case A when both tanks are one-half full (A]_) and when 
the fuel height equals one-half the radius (A 2 ) and for case B when the 
fuel heights in both tanJjs are 2 inches (B 2 ), 3 inches (B^), and 
k- inches (B 4 ). 

In case A the tanks are spherical, somewhat over 4 feet in diameter, 
and centered on the body axis approximately 4 feet in front of and behind 
the airplane center of gravity. The flight conditions are given in 
table I. The fuel weight in the half -full condition is approximately 
25 percent of the total weight. 

In case B the tanks are spherical, centered in the plane of sym- 
metry slightly less than '5 inches below the airplane axis and 4 inches 
in front of and behind the model center of gravity. The diameters are 
8 inches and the total fuel weight in the half -full condition is 
approximately 46 percent of the. total weight. The flight conditions 
are the same as in reference 2. 
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Motions were calculated for certain subcases of the basic cases 
which were obtained by varying certain significant parameters. By 
comparing the resulting motions an attempt was made, to evaluate the 
effect of varying such factors as amount of fuel^ position of tanks, 
and relative natural fuel and airplane frequencies on the distirrbance 
of the airplane motion caused by fuel motion. Also various initial 
conditions were considered to show the effect of initial conditions on 
the resulting motion. In some cases an initial disturbance in sideslip 
was assumed, and in other cases an initial fuel dist\irbance was assumed. 
An initial sideslip of 5° and a fuel displacement of 10° were arbitrarily 
chosen as standard. Since the eqiiations are linear, multiplying the 
initial displacements by a common factor simply multiplies the resulting 
motions by the same factor. In this connection it must be emphasized 
that, if at any time the calculated motion in any degree of freedom 
becomes too large to satisfy the equations, the following motion is 
meaningless For example, if a 5° displacement in p gives rise to a 
fuel motion much greater than 30° to i|0°, the accompanying p motion is 
meaningless because the assumption of small displacements is violated. 
However, if multiplying the fuel motion by some arbitrary factor, for 
example 2/5^ will bring its peaks down to less than 30° to 4o°, then 
the p motion resulting from an initial P disturbance of 2° can be 
obtained by simply multiplying the previous p motion by 2/5 also. The 
effect of large fuel displacements must be discussed qualitatively. 

The motion in sideslip and the motion of the two fuel pendulms in 
the various subcases are shown in figiares 2 to 12. Comments on these 
motions are presented to facilitate interpretation of the figures. The 
period of each fuel pendulum is called the natural fuel period. The 
period and damping of the airplane, the fuel being disregarded, are 
called the natural airplane period and damping. 

Case Aj,- The natural fuel periods for case A^^ (half -full tank) 

are approximately 1.66 seconds and the natural period of the airplane 
alone is 1.40 seconds. Damping to half -amplitude occurs in 2 cycles. 

The motion in figure 2 is for initial Pq = 2°. The early p motion 
seems to have more damping than the natural airplane mode. The disturb- 
ance arising from the fuel modes is evident after 2 cycles. That the 
irregular residual oscillation of amplitude l/4°- to l/2° is essentially 
due to the fuel modes is seen from the fact that the dominant period in 
the later motion is approximately 1.6 seconds. Notice that in this case 
a 5° initial Pq would almost immediately cause fuel displacements of 
over 80°, so that the following motion would be radically changed. 

The motions shown in figure 3 for {ii)o = (^2)0 = are quite 
regular and indicate one dominant mode in each motion. The fuel period 
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is 1.6 seconds. The airplane period starts at 1.4 seconds and builds 
up to 1.7 seconds and averages 1.6 seconds. The amplitude of the side- 
slip motion is very small. The largest such motion which could occur 
for this type of disturbance would be for Initial = (^2)0 ~ 

and would give 3 amplitude slightly more than 0.1°. 

The small amplitude of the sideslip motion in figure 3 was sur- 
prising. It was conjectured that for this fuel configuration the fuel 
displacement (?i)o = 10° (^2)0 “ ”10°; corresponding essentially 

to an initial yawing moment, might be more effective in inducing an 
airplane oscillation. (See fig. 4.) Apparently, this configuration is 
more effective inasmuch as the sideslip motion now builds up to an 
amplitude of approximately 0.4°. The energy necessary to induce this 
considerable "snaking" type of oscillation seems to be obtained Initially 
from the rear-tank motion, \dilch is in the proper phase relation with 
the sideslip motion to feed energy into it at the start of the motion. 

The fuel periods in the regular motion are slightly over 1.6 seconds. 
The airplane period increases from 1.5 to 1.8 seconds and has an average 
period of 1.6 seconds. 

Case A2 .- The natural fuel periods in case A2 are approximately 

1.52 seconds and the natural airplane period is 1.49 seconds. Damping 
to half -amplitude occurs in 1.3 cycles. 

, The motion in figure 5 is for initial pQ = 0.5°. The early side- 
slip motion seems to be of greater damping than the natural airplane mode. 
The residual airplane motion arising from the fuel modes sets in very 
quickly and is a regular unstable motion of very large relative amplitude, 
with a period of. approximately 1.7 seconds. Both fuels start with a 
period of approximately 1.5 seconds, which increases to I.7 seconds. 

In figure 6, = -10° and (^2)0 = 10°. The sideslip builds 

up to a fairly regular oscillation of 0.4° amplitude, with the period 
increasing from 1.5 seconds to 1.7 seconds. The fuel motion has a 
period somewhat under 1.7 seconds, with amplitude quickly building up 
to the limits where splashing must become important. 

Case B2.- The sideslip motion shown for case B2 appears to be a 

normal damped oscillation for the first 4 seconds' (fig. 7, 3^ = 5°), 

but then the peaks show a slight irregularity instead of damping smoothly. 
The period of the early sideslip motion appears to be somewhat over 
0.9 second and the motion damps to half -amplitude in less than 2 cycles. 
This motion is very close to the undisturbed airplane mode ^period of 

0.92 sec and damping constant of 1^ cycles j. The natural fuel, periods 
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are O.6I second. The fuel motion is very irregular and obviously con- 
tains considerable amounts of at least two characteristic modes. 

Cases and Bi|. - The motion in cases B^ and Bj^ shown in figures 8 

and 9 , respectively, for 3^ = 5° is very much the same as in the pre- 
vious one, except that the disturbance of the airplane mode in the 
sideslip motion appears somewhat more pronounced as the amount of fuel 
increases. 

Figure 10 (case B4, (^i)o = (^2)0 = shows that the motion in 

sideslip resulting from the fuel displacement is much more irregular 
than in case A. The dominant mode corresponds to a fuel frequency, but 
apparently' the airplane mode is present with considerable amplitude. 

The maximum oscillations are approximately ±1/4°. The sideslip motion 
in this case was much more irregular than for the corresponding initial 
conditions in case A. It was conjectured that this might be caused by 
the fact that in this model both tanks are below the X-axis, so that the 
coupling of the fuel motion with yawing and rolling motions does not 
have the same phase relationship as in case A where one tank is above 
and one below the X-axis. Therefore, in case B^^^ shown in figure 11, the 

front fuel tank was assumed to be above the X-axis, all other conditions 
remaining as in case Bi^., In this case the general type of motion does 
seem to resemble that in figure 3* The sideslip, which builds up 
to ±1/4°, shows a snaking at the fuel frequency. 

In case A the fuel natural periods are very close to the airplane 
period. In case B, however, the fuel period is approximately two-thirds 
of the airplane period. In case shown in figure 12 for 3 q = 5°, 

the value of Cnp of the model has been arbitrarily changed to give 

the model a period very close to the fuel period of 0.66 second. Com- 
parison with figures 2 and 5 shows that the motion in this case is very 
much like the motion in case A. • 1 

Transverse accelerations .- In evaluating pilots* reactions to 
snaking oscillations, the magnitude of the transverse accelerations 
involved in the oscillations has been found to be an important factor. 
Acceleration amplitudes above 0.025g are found to be bothersome, and 
amplitudes above 0.08g are considered very unsatisfactory. Calculations 
of the transverse accelerations involved in several of the previous 
motions were carried out. The magnitudes of the acceleration peaks in 
the residual oscillations were found to, be approximately 0.04g to 0.05g, 
The actual motions are not shown since all the airplane oscillations are 
essentially of the same type as the 3 motions. 
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DISCUSSION OF RESULTS 


The motions described in the preceding section are sufficient to 
give a fair picture of the types of possible fuel effects. Moreover, 
since each motion is just a superposition of the characteristic modes 
of the total system, these motions are often easier to londerstand if 
the characteristic roots are known. Physically it is clear that the 
characteristic modes will not differ much from the natural (uncoupled) 
modes when the interaction between the airplane and fuel is small. 
Comparison of figures 2; 5^ and 12 with figures 7^ 8^ and 9 indicates 
that the interaction between airplane and fuel is strongest when the 
frequency of the airplane is close to that of the fuel, as might be 
expected from comparison with the resonance phenomena exhibited by an 
oscillator driving a system at its natural frequency. For this reason 
the characteristic modes of the total system were calculated in cases A2 
and where the frequency ratio between airplane and fuel natiiral 

frequencies was practically unity. The. natxiral modes are given for 
purposes of comparison. The results are given as follows in terms of 
periods and times to halve or double amplitudes, in seconds: 

Case A2, natural modes: 

= 1.49, T1/2 = 1.91; Pi = P2 =1.52 

Case total system: 

P^ = 1.47, Tj^/g = ^1 = 1-29, T^/2 = 1-34; P2 = 1.67, T^ = 4.56 

Case natural modes: . 

Pq = 0.66, T1/2 = 1.29; Pi = P2 =0.66 
Case total system: * 

P^ =0.65, T^^2 = 30; Pi = 0.54, T2 = 1.5; P2 = O.71, T1/2 =0.48 


In these cases of large interaction it is difficult to identify 
one of the characteristic modes as the airplane mode. The character- 
istic mode in which the period is changed least from the natural air- 
plane period has been called the airplane characteristic mode.' However, 
in figures 5 and 12 this mode is not obviously the dominant one, as the 
airplane mode is in figures "J , 8, and 9. 
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The most important effects to be noted in these particular cases, 
where the fuel and airplane frequencies are eqxxal, are that a charac- 
teristic mode which is very lightly damped with a frequency close to 
the natural airplane frequency exists and that an unstable mode appears. 
In connection with the first of these effects, it would seem that, 
theoretically, certain initial conditions might be found that would 
excite mainly this lightly damped mode in the characteristic solution 
for 3, so that the resulting motion would be a typical snaking. Of 
course the required initial conditions might or might not be practical 
ones. 

The total characteristic modes were also calculated for cases 
and to investigate the changes in motion caused by a hypothetical 

shift of one of the fuel tanks. The natural modes and the characteristic 
modes of the total system are given as follows for comparison: 

Cases B|^ and natural modes: 

= 0.84, T 1/2 = 1*6^; ^1=^2= 

•N 

Case total system: 

Pa = 0 . 88 , T 1/2 = I. 9 I; Pi = 0 . 63 , T 1/2 = 5.45; P 2 = 0.53, T 1/2 = I .61 

Case total system: 

Pa = 0.86, T1/2 = 1 . 1 ^; Pi = 0.61, T1/2 = 19 . 5 ; P2 = 0.60, T2 = 10.2^ 

It is interesting to note that both fuels in figure 10 and also in 
figure 11 seem to follow the. more stable fuel mode in the part of the 
motion shown. The p motion in figure 11 seems to show the effect of 
the unstable mode. It appears in this case that, when the fuel tanks 
are in front of and behind the center of gravity, the configuration with 
one tank above and one below the X-axis gives rise to an unstable mode, 
whereas the configuration with both tanks below the X-axis makes both 
fuel modes stable. 

A comparison of figures 2 to 4 or figures 5 and 6 clearly shows 
that the initial conditions can have a very important effect, since the 
least stable mode does not necessarily become dominant for a long time. 
The p disturbance was chosen as a typical airplane disturbance. On 
the other hand the disturbance of the fuels as an Initial condition 
would seem to be completely artificial. However, these motions are 
believed to give a rough idea of the residual oscillations caused by 
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fuel motion, at least insofar as magnitude is concerned, since, if the 
fuels were still displaced after the airplane motion had practically 
died out, the remaining motion might he considered' to he the type caused 
hy a fuel displacement. From this point of view figures 4, 6, 10, and 11 
seem to indicate that residual oscillations of the order of magnitude 
- of 1 / 4 ° to 1 / 2 ° might he expected in these cases. Actually, figures 10 
and 11 would show oscillations of the order of 3/4° for 30° fuel 
displacements. ' , 

Because of the lengthiness of the calculations, only cases A^, A 2 , 
and were carried out far enough to show the residual oscillations 
following a 3 disturbance. It is evident that the motion in case B]^^ 

(fig. 12 ) resembles the motion in case A (figs. 2 and 5) much more than . 
it resembles the unmodified case B motion (figs. 7 to 9). The reason 
for the smaller relative fuel motion in case Bi^^ is probably the fact 

that the relative fuel mass is considerably larger than in case A. In 
case A2 (fig. 5) the residual oscillation dominates the motion almost 
immediately. These results show that the importance of the residual 
oscillation depends mainly on the closeness of the natural airplane and 
fuel frequencies, that is, on the parameter which might be called the 
frequency ratio. Inasmuch as the previous discussion of the character- 
istic modes indicated that the frequency ratio was also the most impor- 
tant factor affecting the characteristic modes of the system, the 
frequency ratio generally can be seen to be the most important, factor 
determining the disturbance of the normal airplane motion caused by the ' 
fuel. Moreover, case A indicates that for spherical tanks the fuel 
frequency may easily be of the same order of magnitude as the airplane 
frequency. Reference 4 indicates that the same is true for rectangular 
tanks and for arbitrarily shaped tanks of reasonable dimensions. Thus, 
even though the residual oscillations might occirr at fuel frequencies, 
these frequencies would not be distinguishable from the normal airplane 
frequency in the cases where the fuel effect is most pronounced, since 
in these cases the frequency ratio approaches imity. 

The effects of unstable modes cannot be understood without con- 
sidering the nonlinear effects due to splashing of the fuel. For linear 
systems the presence of an unstable mode would imply that the total system 
is unstable. It has been shown, however, that in an actual motion if the 
coefficient of the \mstable mode in the solution for the airplane motion 
is very small compared with the coefficient of one of the stable modes 
then the unstable mode will not appear in the early part of the motion^. 

Now even in cases where the interaction is weak, one of the fuel modes 
(with no natural damping assumed) may be unstable. In such cases the 
unstable mode in the airplane motions will be relatively very small, 
while the unstable mode may be dominant, in one of the fuel motions. 

Then this fuel motion may become very large before the effect on the 
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airplane motion can be observed. At this point splashing will set in 
and the equations of motion no longer hold. Actually the energy lost 
in splashing may damp out the fuel motion so that the effect of this 
fuel motion - and therefore of the unstable mode - may never appreciably 
appear in the airplane motion. 

In case however, the interaction is large, and figure 5 shows 

that the imstable mode does soon become very important in the 3 motion. 
Here again the effects of large-amplitude fuel motion and, in particular, 
the nonlinear effects of fuel splashing should be considered qualitatively. 
In cases A^_ and A2; because of the strong interaction, a sideslip of 5° 

woxxld very quickly give rise to fuel displacements of the order of magni- 
tude of 90° or greater. Physically, it is clear that the airplane motion 
is feeding energy into the fuel motion in this early part of the motion. 
Because of the large amplitudes of the fuel motion, a large component of 
this motion is in the vertical direction, so that a considerable part 
of the energy in the fuel motion will go into creating longitudinal air- 
plane motion. Because of the symmetry conditions previously cited this 
energy will not be fed back into lateral motion, so that the fuel tends 
to stabilize the lateral motion by feeding some of the energy from the 
lateral into the longitudinal motion. Also the txirbulence due to 
splashing will absorb energy which will then be lost altogether from 
the motion. 

One conclusion which can be drawn from this discussion, is that the 
motion following large disturbances may be more stable than that following 
small disturbances. For example, the Initial distxirbances in figures 2 
and 5 have been adjusted so that during the motion shown there is little 
splashing, and considerable residml oscillations are shown. But if the 
initial distiirbances in these cases had been 5° or more in 3, the 
energy lost in splashing in the early motion would possibly be so great 
that the residual motion in 3 would be smaller than that shown here. 

This might explain why some airplanes which definitely showed trouble- 
some fuel oscillations were reported to be more stable in conditions 
of large atmospheric t\irb\ilence than in slightly turbulent atmospheric 
conditions. 

A more important conclusion is that the fuel can be used to sta- 
bilize the airplane motion by introducing turbulence, by use of appro- 
priate baffles for instance. This can be seen by noticing that the 
early part of the 3 motion in figures 2 and 5 is very stable. But 
if most of the energy fed into the fuel in this part of the motion were 
converted into heat through turbulence, then as has been pointed out 
this energy could not be fed back into the airplane motion and the 
residual oscillations would not appear. Since the amoiint of energy lost 
in turbulence cannot be calculated analytically, it woxild seem that an 
experimental investigation of the effects of honeycomb or other 
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t\arbulence-lnducing baffles on the airplane stability would be desirable, 
especially in cases where the airplane and fuel natural frequencies are 
approximately equal. 

Finally it is possible in a strongly \instable case of resid\ial 
oscillations, such as shown in figure 5 ^ that the fuel may lose just 
enough of its energy in splashing to reduce its amplitude to where the 
motion is again smooth. Then, because of the instability of the system 
for small motions, the amplitude might again begin to build up. In . 
this way continued oscillations of a more or less regular nature would 
occ\ir when the calculations neglecting splashing show unstable motion. 
This result is important because it shows that somewhat irregular small- 
amplitude oscillations can be expected when the ratio of the airplane 
natural frequency to the fuel natural frequency approaches unity, even 
for moderate fuel masses of the order of one-tenth the total mass or 
less. 

EFFECTS OF ASPHERICAL TANKS 


The calculations have been carried out for rigid-body motion in 
spherical tanks only. Actually this assumes that for small oscillations 
the fundamental wave motion in spherical tanks approximates rigid-body 
motion. This approximation only applies when the tanks are one-half full 
or less. This restriction is not too serious, however, since the fuel 
motion will generally have its greatest effect in this range. 

It is important to note that the potential energy of the fuels is 
simply the potential energy of a set of harmonic oscillators located at 
the positions R^.. Thus, the same general analysis will apply whenever 

the fuel motion in the tank can be represented in terms of harmonic 
oscillators with given effective mass and spring constant. Reference 4 
has already been mentioned as obtaining such a representation for the 
fundamental mode of a rectangular tank. Usually the fimdamental mode 
will be the most important and involve the greater effective mass. It 
is conceivable that for long tanks the second mode might be of a fre- 
quency closer to that of the airplane, and in that case might be more 
Important. In such a case each mode might be represented by a separate 
oscillator. As has been pointed out, the damping is mainly due to 
turbulence and will be more important for aspherical tanks. For small 
motions, however, the damping may still be neglected. 

From the general derivation of the equations of motion, the most 
important result was the effect of fuel distribution on the coupling 
of lateral and longitudinal motions. It is plain that these results 
apply strictly only to spherical tanks. Consider for example a tank 
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of triangular plan, form located in the symmetry plane and oriented sym- 
metrically with respect to this plane. Because of the symmetrical 
orientation it can be seen that, although sideslip motion will give 
rise to forward and rearward forces (because the pressure forces are 
normal to the diagonal surfaces), the forward and rearward motion will 
give antisymmetrical lateral forces which will cancel. In this case 
the result would.be to feed energy from the lateral into the longitudinal 
motion; this condition would be favorable since the longitudinal motion 
is generally well-damped. For an vinsymmetrically oriented tank of this 
type, energy could be fed back again from the longitudinal to the 
lateral motion and the problem would be quite complicated. 

In general, the-results on coupling for spherical tanks would be 
valid to first order for such symmetrical plan forms as the rectangular 
or the diamond-shape ones. For any simple symmetrically oriented shape 
in the plane of symmetry, a loss of energy from the lateral to the 
longitudinal motion might occur. This condition would be favorable. 
Finally, for tanks outside the plane of symmetry the same considerations 
would be valid if the tanks were symmetrically placed with respect to 
the symmetry plane and symmetrically shaped with respect to the plane 
through the tank center parallel to the symmetry plane. 


CONCLUSIONS 


The following conclusions may be drawn from the theoretical analysis 
presented: 

1 . Considerable disturbances of the normal airplane motion can be 
caused by fuel motion. 

(a-) The most important factor determining the effect is the 
ratio between fuel and airplane frequencies. When these are equal 
even moderate amounts of fuel (one-tenth the total mass or less) 
may cause considerable disturbances. 

(b) The most usual type of disturbed motion is a somewhat 
irregular small -amplitude oscillation and the type of motion is 
strongly dependent on the initial conditions. 

(c) The effects of splashing will be to make the motion 
more stable, and the loss of energy in fuel turbulence may make 
it possible to increase the stability by artificially intro- 
ducing turbulence in the fuel. 



MCA TN 2280 


35 


2. The fuel motion may cause coupling between lateral and longi- 
tudinal motions. 

3 . The derivation of the equations of motion for spherical tanks 
may be applied to any tanks where the fuel motion may be represented 
in terms of harmonic oscillators. ■ 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., November 21, 1950 
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APPENDIX A 

A STEP-BY-STEP SOLUTION OF THE EQUATIONS OF 
MOTION IN MATRIX NOTATION 


A matrix method for solving the equations of motion is given in 
reference 9* The first step in this method is to reduce equations (19) 
to a set of first-order equations by introducing as new variables D 0 ^ 
Dt, D^ 3 ^, and D^ 2 - S’® described in reference 9 for D0 and D^^. This 

transforms the equations of motion into a set of nine linear first-order 
equations in the nine variables. Iii matrix notation the equations may 
then be written as follows when there are no applied forces Cy^ Cjj, 
or Cj: 


A(Dq) + Bq = 0 (A1) 

where A and B are ninth-order square matrices and q is the column 
matrix (or vector)^ the elements (components) of which are the nine 
variables. In partitioned form^ 


A = 

I 4 IO 45 

B = 

045 I -14 








where is the identity matrix of fourth order; 0 ]^^ and 0 ^]^ are zero 

matrices of order (4x5) and (5 X 4)^ respectively; and 


A55 = 


^»b 


0 


■^‘b2^2/^ 

■t^bi^l/^ 

0 

0 



^^b%Z 

\ 

^^b2^2^Z2/^^ 



l^b%Z 

0 

Kl^ 

Ezi/b 

-Exi/t> 

OJ 

0 


-Exs/b 


-1 
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Equation (Al) may be solved for Dq by multiplying through by the 
inverse of A: 

Dq = -A"^Bq = Pq ' (A2) 

Equation (A2) shows that^ when there are no applied forces^ the dif- 
ferentiation operator with respect to nondimensional time may be 

replaced by premultiplication with the matrix P = -A“4. For example, 

note that D^q = D(Dq) '= P(Dq) = P^q and that similar relations would 
result for higher powers. 

Now, by Taylor's expansion, the value of q( s + As) may be obtained 
from the value of q(s) by the series 

q(s + As) = q(s) + ^ Dq(s) + 


I)2q(s) + . . 
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By use of equation (A2), this equatidn can he written 

q(s + As) = {19 + ^ P + p2 + . . riq( s) (A3) 


which is the fiindamental recurrence relation used in the step-by-step 
calculation. In equation (A3)^ is the ninth-order identity matrix. 
The set of initial disturbances q^ being given, the magnitude of the 
step As will determine the munber of powers necessary to obtain a 
given accuracy in the solution. Because of the relatively high fuel 
frequencies and because it was desired to obtain the motion to a rather 
large nvunber of periods with reasonable accuracy, the series in equa- 
tion (A2) was used to the sixth power with Zis approximately I/20 of 
the airplane period. -Thus, the matrix relation (A3) was 


q(s + As) = Qq(s) 


(A4) 


where 


Q H I9 



(As)^ p2 
2 ’. 


+ 



p6 
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TABLE 'I 

STABILITY DERIVATIVES ATO MASS AITO GEOMETRICAL 
' CHARACTERISTICS FOR TWO CASES CONSIDERED 


Parameters 

Case A 

Case B 



-l.Oi+2 

-0.80 



0.17 

^0.17 

"^3 * • 

-0.126 ■ 

- 0 .l 4 

^np 

-0.01552 

^-o.o4o 

"V 

-0.3^2 

-0.30 

V ' 

-0.28 

-0.16 

C^r * 

0.0796 

0.30 

Weight of airplane alone, lb .... 

6970 

11.25 

S, sq ft 

130 

2.67 

b, ft 

28 

4.0 

7 o> deg ..... • . ; . . 

0 

-11 

p, slugs/cu ft . . 

0.00136 

.6.002378 

Uq, ft/sec 

704 

61.5 to 74.25 


®In case B]^^, 0 ^^ = 0.29. 

^Actually, slightly different values of C^p were used for each 
of the subcases of case B. 












MCA TN 2280 


4 l 


TABLE II 

FUEL DEPENDENT PARflMETEES 


Parameters 

Al 

A2 . 

B2. 

®3 

B4 

Radius of forward tank, ft . • 

2.15 

2.15 

0.333 

0.333 

0.333 

Radius of rear tank, ft , * . 

2.12 

2.12 

0.333 

0.333 

0.333 

Fuel height in forward tank, 
ft 

2.15 

1.075 

0.1667 

0.250 

0.333 

Fuel height in rear tank, ft . 

2.12 

1.06 

0.1667 

0.250 

.0.333 



3.5 

3.5 

0.333 

0.333 

0.333 

1^X2 > 

- 4.1 

- 4.1 

-0.333 

-0.333 

-0.333 

ft 

0.0123 

0.0123 

0 . 4 o 7 

0 . 1 lO 7 i 

0.407 


- 0.0144 

- 0.0144 

0.407 

0.407 

0.407 

Weight of forward fuel, lb . . ■ 

l 48 o 

462 

1.51 

3.07 

4.85 

Weight of rear fuel, lb ... 

1068 

334 

1.51 

3.07 

4.85 

Total weight of airplane and 
fuel, lb 

9518 

7766 

14.27 

17.39 

20.95 

Cl . . 

0.19 

0.19 

1.062 

1.178 

1.174 

Ix^ slug-ft^ 

1360 

1360 

0.1081 

0.1243 

0.1428 

I2, slug-ft^ 

7340 

7708 

0.2098 

0.2206 

0.233 

Ixz, slug-ft^ 

262 

277 

0 

0 

0 



0.806 

1.45 

0.225 

0.174 

0.125 

^ 2 ^ ft 

0.795 

1.43 

0.225 

0.174 

0.125 

Ij^, slug-ft^ 

85.2 

4 o.o 

0.00315 

0.00507 

0.00668 

I2 , slug-ft^ 

59.7 

28.1 

0.00315 

0.00507 

0.00668 
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;ure 3.- Present-day high-speed airplane with tanks one-half full 
(case A^). Sideslip and fuel motions following initial fuel 
disturbance, (^1) = =10°. 
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Figure 6.- Present-^y high-speed airplane with fuel heights one-half 
the tank radius (case A2). Sideslip and fuel motions- following 
initial fuel disturbance, (?i)n = -10°, ( ^„") = 10°. 
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Figure 7.- Free-flying, airplane model with fuel heights one-half the 
tank radius (case B2). Sideslip and fuel motions following initial 
sideslip, 3o = 








.eslip and fuel motions following initial fuel 



Figure 11.- Free- flying aiiplane model with tanks one-half full and 
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